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We study conductance through a quantum dot under Coulomb blockade conditions in the presence 
of an external periodic perturbation. The stationary state is determined by the balance between the 
heating of the dot electrons by the perturbation and cooling. We analyze two cooling mechanisms: 
electron exchange with the cold contacts and emission of phonons. Together with the usual linear 
Ohmic heating of the dot electrons we consider possible effects of dynamic localization. The com- 
bination of the abovementioned factors may result in a drastic change of the shape of the Coulomb 
blockade peak with respect to the usual equilibrium one. 
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I. INTRODUCTION 

At low temperatures electronic conduction through a 
quantum dot weakly coupled to the contacts is governed 
by the Coulomb blockade effect pj - suppression of trans- 
port due to the energy cost of changing the number of 
electrons in the dot. Efficient conduction through such 
a dot is possible only when the electrostatic potential of 
the dot, controlled by external gates, is tuned to a spe- 
cial value where the Coulomb energies of the states with 
N and N + 1 electrons in the dot are close for some N. 
As a result, the linear response conductance exhibits a 
sharp peak as a function of the gate voltages. Theory of 
the Coulomb blockade in equilibrium is well developed 
by now [2(. 

In the last few years several experiments have been 
done on quantum dots under an external ac perturba- 
tion Q . Under these non-equilibrium conditions the elec- 
tronic temperature of the dot is no longer determined by 
the external cryostat, but by the balance between heat- 
ing by the ac perturbation and cooling due to various 
mechanisms. At sufficiently low temperatures cooling is 
dominated by simple electronic exchange between the dot 
and the contacts (the latter are assumed to be maintained 
at a constant low temperature determined by the cryo- 
stat). In this case, as the gate voltage is tuned away from 
the Coulomb blockade peak, the cooling rate changes, so 
does the electronic temperature, thus changing the peak 
shape with respect to the equilibrium one. This simple 
qualitative consideration poses the problem, which is go- 
ing to be studied in detail in the present work. 

Another motivation to study these effects is the search 
for experimental signatures of dynamic localization (DL). 
Experimental observation of DL in trapped ultracold 
atoms in the field of a modulated laser standing wave Q| 
provided a solid ground for the preceding extensive the- 
oretical studies of the kicked quantum rotor 00. In 
a recent publication 0] we have shown that an analo- 
gous suppression of the energy absorption is possible for 
a solid-state system - a chaotic quantum dot under an ac 
excitation, e. g. like those used in experiments of Ref . 0] , 



which makes the question about the possibility of obser- 
vation of DL in a quantum dot highly relevant. If one 
wishes to detect this effect by transport measurements, 
the Coulomb blockade regime is the most suitable, since 
it is in this regime that the transport is sensitive to the 
internal state of the dot, while for an open dot, when 
electron-electron interaction can be neglected, the con- 
ductance is insensitive to the electron energy distribution 
in the dot 00. 

First, consider the standard picture of heating by an 
ac perturbation. Let the single-electron mean level spac- 
ing 5 in the dot be small enough. Then, if an exter- 
nal time-dependent periodic perturbation with the fre- 
quency u) is applied, the total electronic energy E in the 
dot (counted from that of the ground state) grows lin- 
early with time as described by the Fermi Golden Rule: 
E{t) = TujH/S = W t. The probability of each single- 
electron transition per unit time, denoted by T, measures 
the strength of the perturbation ^(| . The criterion of va- 
lidity of the Fermi Golden Rule is S <C T, and T <C w is 
also assumed (H = 1). This picture corresponds also to 
the classical Ohmic absorption by a small particle made 
of a metal with large conductivity a u). 

This picture (hereafter referred to as Ohmic absorp- 
tion) is valid provided that each act of photon absorption 
by an electron is independent of the previous ones; how- 
ever, for a discrete energy spectrum this turns out not to 
be the case. After many transitions the absorption rate 
decreases due to accumulation of the quantum interfer- 
ence correction 0, so that after a time i* ~ T/S 2 the ab- 
sorption is completely suppressed. This effect was named 
the dynamic localization in energy space; the effective 
electronic temperature (the characteristic spread of the 
electron distribution function), reached by the time t*, 
T* ~ Toj/8, plays the role of the localization length. Note 
that DL has nothing to do with the saturation of absorp- 
tion by a pumped two-level system, as in our case the 
spectrum is unbounded. DL is the consequence of level 
discreteness: at S — > it takes longer time for the DL to 
develop, and for the continuous spectrum there is no DL. 
Since this effect drastically modifies the heating rate, the 
stationary state of the dot is strongly affected. 
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The considerations of Ref. pj were based on random 
matrix theory description of the single-particle properties 
of the dot. This description is valid provided that all 
energy scales in the problem are small compared to the 
Thouless energy Exh (defined by the order of magnitude 
as the inverse of the time required for an electron to travel 
across the dot and thus to randomize its motion due to 
scattering off the dot boundaries). For the dot to be in 
the Coulomb blockade regime, the effective temperature 
should be also smaller than the dot Coulomb charging 
energy E c . Thus, in the following, the hierarchy of scales 
J « T « w « T, « E T h, E c is assumed. Note that in 
the random matrix theory one can neglect multiphoton 
processes as they are of the order of the inverse matrix 
size. 

Possible cooling mechanisms for electrons in the dot 
are (i) electron exchange with the contacts, and (ii) en- 
ergy exchange with the phonon subsystem. Both elec- 
trons in the contacts and phonons in the dot are as- 
sumed to be maintained at a constant temperature T 
determined by the cryostat. In the following we analyze 
the interplay of the abovementioned effects in heating 
and cooling, and see how they affect the shape of the 
Coulomb blockade peak. In a short preliminary version 
of this study we have considered only the first cooling 
mechanism . Here we include cooling by phonon emis- 
sion which, to the best of our knowledge, has been little 
studied for a quantum dot. 

The paper is organized as follows. In Sec. II we analyze 
the heating and discuss how it is affected by dynamic lo- 
calization. Sections III and IV are dedicated to a detailed 
analysis of the two cooling mechanisms. In Sec. V we 
consider the resulting stationary state and the Coulomb 
blockade peak shape. Finally, in Sec. VI we summarize 
the main results. 



II. HEATING BY AC PERTURBATION 

In the Ohmic regime the energy absorption by elec- 
trons is linear in the field intensity and given by Wq = 
Tuj 2 /5 (we remind the reader that T is a measure of the 
microwave field intensity, equal to the probability per 
unit time of a single one- photon transition). The same 
expression can be obtained from simple classical argu- 
ments considering a small particle made of a metal with 
a large finite conductivity a 3> co- 
in the regime of the strong dynamic localization the 
absorption is no longer given by the simple Ohmic ex- 
pression. For non-interacting electrons in a closed dot 
the absorption becomes completely suppressed by inter- 
ference corrections that develop in a characteristic time 
i* ~ r/<5 2 , and the effective temperature of the elec- 
trons, reached by that time, is T* ~ Tui/6. In the 
presence of weak dephasing processes with the dephas- 
ing rate 7^ <C there is a residual absorption with 



the rate given by 

W in ~ W ol4 ,U —T^^j-. (1) 


If the dephasing is too strong, 7^ > the dynamic 

localization is destroyed and W- ln ~ Wq. 

The expression was justified in Ref. ^2 f° r the 
dephasing due to electron-electron collisions. The main 
condition of its applicability is that dephasing should be a 
sequence of distinct phase-destroying events with average 
frequency 7^,, rather than phase diffusion, in which case 
the dephasing rate roughly coincides with the quasipar- 
ticle relaxation rate: 7^, ~ j qp . This is certainly correct 
for the case of electron escape to the contacts, since in 
this case the electron is effectively replaced by another 
one with an absolutely random phase. This is also true 
for electron-electron and electron-phonon collisions, since 
the typical energy transfer during a collision is of the or- 
der of the (effective) electronic temperature in the dot, 
which is large: T 3> "f qp (this inequality follows from 
T* > l/U due to 6, and from T > T*). 

Once the condition 7^ ~ j qp is verified, the follow- 
ing consideration can be applied. As the collisions are 
rare (7 gp i* -C 1), the electrons spend most of the time 
in the states localized in energy space, having definite 
phase relationships. When at some moment the phase of 
some electron is destroyed, its wave packet starts spread- 
ing along the energy axis. It localizes again after the 
time ~ t*, in the meantime spreading by ~ T*. Thus, 
the ac driven dynamics following the collision leads to 
a change of the total electronic energy of ~ T* per col- 
lision. The sign of this change is, however, arbitrary, 
because a periodic perturbation can equally cause tran- 
sitions up and down the spectrum. Only the presence 
of the filled Fermi sea below (i. e., an energy gradient 
of the electronic distribution function) makes absorption 
the preferred direction, which means that if the electronic 
temperature T 3> T», the energy absorbed per collision 
is on the average ~ T^/T rather than T». The effective 
number of electrons that can participate in a collision is 
~ T/S (due to the degenerate Fermi statistics). During 
the time interval ~ 1 /j qp each of these electrons partic- 
ipates in one collision, so the total number of collisions 
per unit time is ~ (T/S)j qp . This gives the energy ab- 
sorption rate W- m ~ (T 2 /T)(T/S)^ qp , which is exactly 
Eq. 0. 

The same can be seen from an alternative argument. 
After each collision the electron spends the time ~ t* 
absorbing the energy from the microwave field, then it 
stops to absorb (the dynamic localization occurs) and 
waits for the next event (provided that i* <C l/7</>). Thus 
the absorption rate of the whole system is given by the 
simple weighted average: W m ~ WoJ^t*, which is again 
Eq. 0. 

An important point is that dephasing rate, generally 
speaking, depends on the electronic temperature, which 
results in a temperature-dependent absorption rate in 
the DL regime (Fig. 0. The temperature, in turn, de- 
termined by the balance between energy absorption and 
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T* T, a.u. 



FIG. 1: A schematic view of the dependence of the absorp- 
tion rate W ln on the effective electronic temperature in the 
dot: in the dynamic localization regime the absorption is due 
to dephasing, so it is temperature-dependent; when the tem- 
perature becomes high enough, the dephasing destroys the 
dynamic localization and the absorption is Ohmic. 

cooling. This feedback leads to a non-trivial dependence 
of the characteristics of the stationary state on the con- 
trol parameters, which will manifest itself in a change of 
the Coulomb blockade peak shape, as will be shown be- 
low. The absorption itself becomes nonlinear with the 
field intensity through the dependence of 7^ on Wq- 



If all gates have the same voltage, then (up to a constant) 
U is given by this voltage. Generally, we will call U the 
reduced gate voltage] it is a natural control parameter for 
the system. If the expression in the brackets is of the 
order of unity and the temperature T <C E c , the conduc- 
tance through the dot is suppressed due to the Coulomb 
blockade. If the gate voltages are tuned so that the ex- 
pression in the brackets is small for some particular N, 
the dot conductance G{U) exibits a peak for these values 
of U. The width of the peak AU ~ T, which can be used 
to measure the temperature of the system. 



A. Sequential tunneling 

When U is tuned to the peak, the main contribution to 
the conductance comes from the leading order of the per- 
turbation theory in the dot-contact coupling. For charac- 
teristic temperatures T> S one can describe the system 
by rate equations of Kulik and Shekhter |l3( . We consider 
these equations for the case when the electron energy dis- 
tribution function in the dot f e is non-equilibrium. Let 
the distribution in the ath contact be Assuming 
the dot to have either N or N + 1 electrons with the 
probabilities pn, Pn+i to have N or N + I electrons on 
the dot (all others are neglected, so p^ +pn+i = 1), we 
can write the rate equation as 



III. COOLING DUE TO ELECTRON ESCAPE 

We characterize the coupling of the dot to the two 
contacts by single-particle escape rates 71 and 72 . When 
they are much smaller than the mean single-particle level 
spacing 8 in the dot, the fluctuations of the total charge 
on the dot are small. If the dot is coupled to several gates 
through capacitances C% and voltages Vi are applied to 
the gates, the electrostatic energy of the dot with N elec- 
trons on it is given by 



E(N) 



e 2 N 2 
2C 



C = J2d, (2) 



where e 2 /(2C) = E c is the charging energy. The energy 
cost of adding an electron is 



U = E(N + 1) - E(N) = 



C 



N 



E 



CM 



(3) 



^ = w£7«//.(i-/S)f- 

Q — 1,2 

-2 PN 7«/(l-/ e )/Sf , (4) 



where the factor of two comes from the spin degener- 
acy. The distributions in the contacts are assumed to be 
Fermi-Dirac ones with the temperature Tq: 



To _ 



1 



e e/T + I 



(5) 



As usual, we require pn and pn+i to be stationary. Shift- 
ing the distribution in one of the contacts by an infinites- 
imal voltage, one obtains the linear response conduc- 
tance G: 



G{U) = — 

7i 



2e 2 7172 Fl(U) F^{U) - F( n (U) F Q 2 ut (C/) + F in (U) F out (U) 



■ 72 



F ta (U)= I (l-f e )f?° u( le, 



[F in (U) 4 
F out (U) 



m 2 

/e(l- 



(6) 
(7) 
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In the equilibrium case, when / e = fj° as well, the last 
fraction in the right-hand side of Eq. |JB} reduces to the 
familiar expression (1/2)(U/T)/ smh(U/T). 

Tunneling events lead to the change in the distribution 
function in the dot. The kinetic equation describing this 
process can be obtained straightforwardly from equations 
of Ref. Hal and reads as 



-QjT = (71+72) X 

(1 - MfT+u F out(U) - W - fllu) FUU) 
F in {U) + F out (U) 

If we introduce the functions 



£ in (U) = Hl-fJfSuede, 



£out(U)= I f e (l-f?°)ede 



(8) 

(9) 
(10) 



and denote by 7 = 71 +72 the total single-electron broad- 
ening, the cooling rate for the dot electrons (the total 
energy loss per unit time) can be written as: 



Want (17) 



7 £ou_ t (U) F in (U) 
S 



£in(U) F out (U) 



F in {U) + F out (U) 



(11) 



From the kinetic equation (jHJ) one can also extract the 
single-particle escape rate for a particle with the energy e: 



lesc = 7(1 - fllu) 



F in (U) 



F ia {U) + F out (U) 



(12) 



In the following we will use the expression for "f esc at 
e = as an estimate. We will also use the Fermi-Dirac 
form for the electronic distribution function with some 
temperature T. This is true only if electron-electron col- 
lisions restore the Fermi-Dirac shape much faster than it 
is modified by other processes. If this is not the case, 
T still gives the characteristic width of the distribution 
function. It is determined by the balance between heat- 
ing by the ac field and cooling considered in the previous 
section. 

We also assume the electronic temperature in the dot 
to be much higher than the temperature of the contacts 



J 



G(U) 

dfe 



4e 7172 



7 



irU 2 5 2 
2 



(the latter can be made as low as ~ 10 mK [Tj|), which 
is true if the pumping power is high enough. Then we 
can set the temperature of the contacts to be zero, which 
allows an explicit calculation in Eqs. ©~dJ [we denote 
x = U/(2T),G = G(U = 0)}: 



-2x1 



F in (U) = T\n [1- 
F ont {U) =F ia (-U), 
G(U) ^ xtanhx 
Go In (2 cosh x) ' 

00 

£ in (U) = -2T 2 J (l-tanhy)ydy, 
£out{U) = -& m (-U), 

W out (U) _ IT 2 2 

(l/S)T 2 12 X + 

X 

2x 

ln(2 cosh a:) 



y tanh y dy , 



lesc{U) 1 
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2 21n(2coshx) 



(13) 
(14) 

(15) 



(16) 
(17) 

(18) 
(19) 



B. Inelastic cotunneling 

At large U S> T sequential tunneling becomes sup- 
pressed exponentially. In this situation both the con- 
duction and cooling become dominated by cotunneling 
- a second-order process whose probability contains an 
additional small factor 7 5 /U 2 . Obviously, only inelastic 
cotunneling |l5j| can contribute to cooling. Elastic cotun- 
neling 16], which does not change the electronic state 
of the dot, contributes to conduction at temperatures 
T < \JE C 8. We will be interested in higher temperatures 
and do not consider this contribution. 

A straightforward generalization of the considerations 
of Ref. to the non-equilibrium case leads to the fol- 
lowing expression for the conductance in terms of the 
electronic distribution functions and the kinetic equation 
for the distribution in the dot: 



/ e _n(l-/e)(l-/J°) 



Of, 



To " 
e'+Q 



de' 



de de' dfl . 



dt irll 2 5 



(1 - / e )/e-n(l - - Ml ~ / £ -a)/J°(l - f? +n ) 



de' dn 



(20) 
(21) 



For a Fermi-Dirac distribution, f e — / e T the integrals can be calculated explicitly for any temperatures T, Tq: 



G(U) 



27re 2 7i72 , 2 , ^2 



3U 2 S 2 



(22) 
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w out (u) = 2 ^,gT (r 4 -^ 4 )- (23) 

The electron escape rate at e = can be extracted from 
the kinetic equation (|21|l : 



7e 



7T 7 



(24) 



If we set, as before, To = 0, we obtain the following 
explicit expressions: 

W ont 7T 3 j/S G 7T 7/(5 7 esc 7T j/5 



(7/<5)T 2 30 x 2 ' G 6 x 2 ' 7 24 x 2 ' 

(25) 



C. Photon-assisted tunneling 

So far the only effect of the AC perturbation we were 
interested in was to cause transitions between single- 
particle states in the dot. The perturbation, however, 
may possess a component Ucoso->£, proportional to the 
unit matrix in the dot single-particle Hilbert space. In a 
closed dot this component does not cause any transitions 
and can be gauged out completely, so it does not affect 
any observables, either single-particle or many-particle 
ones (in particular, it does not affect electron-electron 
collisions). 

However, when the dot is connected to contacts, this is 
no longer the case, as the diagonal component is respon- 
sible for the photon-assisted tunneling 0. This effect 
can be taken into account by replacing the dot electron 
distribution function f € in the above formulas by 



(26) 



where J n is the Bessel function. Photon-assisted tunnel- 
ing will not be important for our considerations if the 
smearing of the distribution function given by Eq. I|26l) 
is much smaller than the thermal smearing. Using the 
asymptotic expansion of J n (z) at large n, this condition 
can be written as 



J n {z) 



max{V»<T. (27) 



The condition uj -C T is authomatically fulfilled if T ^S> S 
and T > T* = Tu/5. As for the condition V < T we note 
that within the N x N random matrix approximation, 
adopted in Ref. [lj, we have (V 2 ) = (1/N)T6 so that 
V -> as N -> 00. 

Besides the random component with zero mean in- 
cluded in the random-matrix treatment, V can have a 
deterministic part. It is given by the spatial average of 
the perturbation potential over the dot volume, and en- 
ters our model as an independent parameter. Thus in 
order to fulfill the condition FcTa special experimen- 
tal care should be taken. 



IV. COOLING DUE TO PHONON EMISSION 
A. General expressions 

Another important mechanism of electronic energy re- 
laxation is emission of phonons. For mesoscopic metallic 
rings with diffusive electronic motion this problem was 
addressed in Ref. 0|. For quantum dots energy relax- 
ation at frequencies smaller than the mean level spacing 
has been considered ; here we are interested in 

the opposite limiting case, 5 being the smallest energy 
scale. Below we estimate the corresponding cooling rate 
for clean (ballistic) quantum dots made out of 2D elec- 
tron gas (2DEG) in a GaAs/AlGaAs heterostructure Q 
and bulk 3D phonons. 

For ballistic dots (whose size L is smaller than the 
elastic mean free path £) one does not need to take into 
account phonon-induced impurity displacements |2l|, so 
the phonon-induced potential felt by the electrons can be 
written in the form 



V(r) = 



rf 3 q 
(2tt) 3 



U(q)e 1 ' 



q,A 



h. C. 



(28) 

where bq \ is the annihilation operator for a phonon 
mode A with the wave vector q. The detailed form of 
the coupling depends on the specific coupling mech- 
anism to be specified below. 

The probability of the electronic transition from an ini- 
tial single-particle state s with the energy e s and the wave 
function ip a (r) to the final state s' with the energy e s > and 
the wave function ip s >(r), accompanied by absorption or 
emission of one phonon, is given by the Fermi Golden 
Rule: 

2 



.abs(om) f ^(r)e ±lqr 



N, 



EsT^a), (29) 



where N^x is the phonon occupation number before the 
transition, and ujq.\ is the phonon frequency; the upper 
sign corresponds to the phonon absorption, the lower one 
- to emission. Introducing the transition rate 

4^) = ^EK-'+^)^-e,)5(e'-£ 8 '), (30) 



averaged over the random dot realizations, we can write 
the kinetic equation for the electronic distribution func- 
tion f e : 

Me',e)(l-f £ )f e ,-w(e,e')f e (l-f £ ,)}^-. 

(31) 

The average rate l|30|) is determined by the electronic 
wave function correlations in the dot: 



n e ,e'(r,r') 



> s (r)V s V)Vv(r')e(r)x 
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x5(e- e s )S(e' - e s >), 



(32) 



averaged over the dot realizations. Then we can write 
the average transition rate as 

w(e,e') = 2^ 2 ^n^(q,q)|Vq,A| 2 x 

x [7V q , A 5(e' - e - w q , A ) + 

+ (JV q , A + l)(J(e'-e + a; q , A )], (33) 

with the Fourier transform defined as 

n e , e ,(q,q') = J n £!£ -(r,r')e^ lqr+4q ' r '(i 3 rd 3 r'. (34) 

Statistical properties of ballistic dots have been ex- 
tensively studied (in Refs. |22Tl23T ]. for a review see 
Refs. For |e - e'| smaller than the Thouless en- 

ergy Exh one can use the following estimate: 



n e , e <(q,q) 



1 



E Th 5 



iin{l, (<7||L) 2 }, 



(35) 



where the factor {q\\L) 2 appears when q\\L <§C 1 (qy is the 
component of the wave vector parallel to the plane of the 
2DEG). As a result, the transition rate w(e, e') depends 
only on the transferred energy u> = e — e' . 

We assume the electronic temperature (determined by 
the balance between heating and cooling) to be much 
higher than the lattice temperature (determined by the 
external cryostat). In this case one can neglect any 
phonon population present, = 0, so only emission of 
phonons can occur, and w(lo) oc 9(lo). For a power-law 
dependence, w(lo) oc uj°'9{lj) 1 and Fermi-Dirac electron 
distribution in the dot © the cooling rate is given by 

Wout= / j ™iow(u)f e (l-f e _ u )oc-^-. (36) 

Obviously, such a power-law dependence can be 
parametrized by a single parameter T p h and written as 
W out = T Q+3 /T p a h +1 . From the kinetic equation (5TJ one 
can also extract the single-particle relaxation rate 7 P h- 
For electrons with the typical energy e ~ T and w(lo) oc 
uj a 9(uj) we obtain 7 ph (T) - 5 (T/T ph ) a+1 . 



B. Specific mechanisms 

To consider specific electron-phonon coupling mecha- 
nisms, we describe phonons in terms of the lattice dis- 
placement operator for each normal phonon mode A: 



1 



2Vp A W q ,A 



-zqr 



(37) 

The displacement of each mode is directed along the unit 
vector e q a- To each mode corresponds some mass which 
is the total mass of the unit cell for acoustic phonons 



or the reduced mass for optical phonons; dividing it by 
the unit cell volume one obtains the corresponding den- 
sity p\. Finally, V is the 3D quantization volume. At low 
temperatures we are interested in, only acoustic phonons 
can be emitted. We approximate their dispersion by 
Lo q = v s q, with v s being the sound velocity, while the 
density p\ coincides with the density of the crystal po- 

Dej ovulational coupling to the acoustic phonons is due 
to the local change of the electronic energy bands under 
strain: 



V« ct (q)=Z 3l iqMl) 



(38) 



where S-,-; is the deformational coupling tensor. In a bulk 
crystal it, generally speaking, depends on the electronic 
wave vector k |2J|. In doped GaAs, when typical elec- 
tronic wave vectors are close to the Brillouin zone cen- 
ter and one approximates the periodic part of the Bloch 
function by that for k = 0, this dependence vanishes and 
Sj/ = ZSji. The leading anisotropic (i. e., dependent on 
the direction of k) correction at small finite k should be 
smaller by a factor of (fca) 2 , where a is the lattice con- 
stant (the first-order in k correction should vanish due 
to the time-reversal symmetry). Hence we can estimate 
its magnitude as ~ (fci?a) 2 S 
electronic Fermi wave vector, and 
density. 

The isotropic (independent of the direction of k) par t 
of the deformation potential is subject to screening [25j. 
The electrons inside the dot can screen the fields with 
wave vectors q\\ down to ~ 1/L. We assume that the 
Fourier components with q« <C 1/L are also screened, 
either by the 2DEG outside the dot, or by the metallic 
gate. Thus, by the order of magnitude, we can use the 
expression for the static (due to v s <C vp) screening by 
an infinite 2DEG, which results in the renormalization: 



na H, where &f is the 
is the 2D electron 



l + l/(?||O a ) 



(39) 



where a s is the 2D screening length (equal to half the elec- 
tronic Bohr radius), and we consider q\\a s -C 1. Thus, the 
effective deformation potential is suppressed by a small 
factor: either by q\\a s for the isotropic part of the poten- 
tial, or by no? for the anisotropic part. 

Piezoelectric coupling to acoustic phonons is due to 
the longitudinal electric field induced by the strain. We 
express the potential in terms of the electromechanical 
tensor e™, which relates the induced polarization to the 
strain tensor: 



V pl ' 



(q) = ^--^«j(q). 



(40) 



Here e is the background dielectric constant of the ma- 
terial. The in-plane piezoelectric field is also subject to 
screening, which brings a factor of q\\ a s . 

The component of the piezoelectric field perpendicular 
to the dot plane is not screened by the electrons. In- 
stead, it affects the confinement and shifts the subbands, 
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which can be viewed as Stark effect. If q\\ ^ 0, the shift 
of the subband depends on the in-plane coordinate and 
represents an additional effective potential felt by the 
electrons. If we assume the confinement of the electrons 
by an asymmetric triangular potential well formed by 
the constant force F on one side and a hard wall on the 
other, the confinement energy e z ~ (h 2 F 2 /to) 1 / 3 , while 
dcz/dF = a z is of the order of the extent of the confined 
state in the z direction. Thus, we can estimate 



l/ btark (q) ~ q z a z VP lezo (q) . 



(41) 



STATIONARY STATE 



A. Ohmic absorption 



First, consider the case of the simple Ohmic absorption 
with cooling only due to the contacts in the sequential 
tunneling regime with the rate given by Eq. (|18|l . At 
small detunings (U <C T) we have 



Wout(T) = J 



n 2 T 2 
12 



U 2 



(44) 



This effective in-plane potential is also subject to screen- 
ing, which brings an additional factor of q\\a s - 
As a result, we can generally write 



V a ~ A 



Vp v s 



(42) 



with A given by 



qa s ^, na 2 E, 



47ree om a„ 



47ree om a s 



, q z a z 



for the screened isotropic deformation potential, 
anisotropic deformation potential, screened in-plane 
piezoelectric field, and the perpendicular piezoelectric 
field, respectively. 

Let us estimate the relative importance of these mech- 
anisms, using the numbers for GaAs from Ref. j2(|. The 
bare deformation potential 5 ~ 10 eV, the screening 
length a s w 50 A, the lattice constant a ~ 5 A. We 
will be interested in temperatures T ~ 0.1-1 K, so we 
indeed have qa s <C 1, and we are in the regime > 1 
(for L ~ 1 fim). For v s q — 1 K we have q ps 3 • 10~ 3 A -1 
(v s 5 • 10 5 cm/s), so qa s ps 0.15. For n = 10 12 cm -2 
na 2 ps 2.5 • 10~ 3 , so the screened isotropic part is more 
important than the anisotropic one. The only indepen- 
dent component of the electromechanical tensor in GaAs 
ef™ w 1.4 ■ 10 7 V/cm, the dielectric constant s ps 13, 
so Aireef£a s / e ps 7 eV, which is of the order of the un- 
screened deformation potential. For the screened poten- 
tial due to the perpendicular piezoelectric field, as typ- 
ically a z ~ 100 A [23, we have a smallness of q z a z . In 
conclusion, contrary to the estimates of Ref. |20j , we ob- 
tain that the in-plane piezoelectric coupling is more im- 
portant than the deformational one. 

As a result, we arrive at the estimate 



so that the stationary temperature is given by 



{l/5)U 2 
8W 



OiU A ) 



(45) 



The temperature T(U = 0) determines the curvature of 
the Coulomb blockade peak at U — 0: from Eq. i|15fl we 
have 



G(U) 
Go 



= 1 - 



1 



U 1 



41n2 T 2 {U = 0) 



(46) 



At large detunings 3> T we can approximate the right- 

2 1 x and write 



hand side of Eq. l(T5|l by \x\e 
W 



{l/5)T 2 



-2N 



T 



U 



ln[( 1 /6)U 2 /(2W )} ' 



(47) 



with the logarithmic precision. It is correct if the loga- 
rithm in the denominator is large, or U ^ T(U = 0). 
This result means that the tails of the Coulomb blockade 
peak have the form: 



2Wn 



G(U) 

Go ~ (j/5)U 2 



In 



(l/S)U 2 
2W 



(48) 



The weak power-law fall-off of the tails is drastically dif- 
ferent from the exponential one occurring in equilibrium: 
G(U)/G = (U/T)/sinh(U/T) 0. The reason for this 
difference is very simple: as the gate voltage is tuned 
away from the degeneracy point, the exchange of elec- 
trons between the dot and the contacts becomes weaker, 
so the cooling rate decreases leading to an increase in the 
temperature and hence in the conductance. 

At large enough detunings the cooling becomes domi- 
nated by the inelastic cotunneling, Eq. I|25|l . rather then 
sequential tunneling, Eq. I|18[). In this regime the dot 
temperature and the conductance are given by 



W out (T) 



(4nee cm a s /e) 2 
pavlE Th 5 



2^6 

7pr . ( 43 ) 



For GaAs p v 5 s w (0.074 eVY* (the density p 
5.3 g/cm 3 ), for a typical dot :3j 5 ~ 1 peV, E Th 
100 fxeV, so we obtain T ph - 0.1 meV - IK. 



T = 



15 W U 2 



G{U) 



2tt 3 ( 7 /5) 2 J ' Go V 3tt 2 U 
The switching to the inelastic cotunneling occurs at 

G(U) j/S 



U ■ 



Wo S -ln 25 -, 

1 1 



Go In (<5/7) 



(50) 
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The logarithmic precision of these estimates, however, 
makes them applicable only for extremely small 7 [such 
that ln((5/7) 3> 1]. In reality if one takes directly the 
expressions (|18f) and H25|) for the cooling rate, for Wo = 
30 fieV 2 w 46 fieV/s, j/8 = 0.2 28] the contribution of 
the inelastic cotunneling starts to affect the stationary 
electronic temperature noticeably [as compared to the 
precision of Eq. JSJl] only as far as V > 1 meV. At 
V = 1 meV the conductance G(l meV)/Go « 0.01, and 
about 18% of it is still due to the sequential tunneling. 

If one takes now into account cooling by phonons with 
the rate l|43|) ■ it sets the upper limit for the electronic 
temperature: T max = (WqT^) 1 / 6 . If the pumping is 
strong enough (or the dot is closed enough), T max <g; 
y/WoSpy, the phonon cooling mechanism dominates, the 
dot temperature is constant and equal to T max for all U, 
so that the Coulomb blockade peak shape is given ex- 
plicitly by Eq. (JTHJl, and its tails - by Eq. (|2*3|l . In the 
opposite limiting case, T max 3> y/WoS/j the electronic 
temperature in the peak region is determined by elec- 
tron escape, and only in the peak tails, when the dot 
effectively becomes more and more closed, phonon emis- 
sion starts to dominate. This will manifest itself as a 
crossover from the 1/U tail l|49|l to the 1/U 2 one given 
by Eq. I|25|) at fixed T = T max . This crossover occurs at 
U ~ h/S)Tl ax /Vm, G/Go ~ (S/lWo/T^ P. We 
plot the tails of the Coulomb blockade peak G(U)/Gq in 
Fig. |21 for three cases when (a) only sequential tunnel- 
ing is taken into account, (b) cotunneling is added, and 
(c) cooling by phonons is present as well [28j. In Fig. [3] we 
plot the electronic temperature T(U) for the same three 
cases. 




100 200 
U, neV 



300 



FIG. 2: Normalized conductance versus reduced gate volt- 
age (Coulomb blockade tail): (a) only sequential tunneling is 
taken into account, (b) cotunneling is added, and (c) cooling 
by phonons is present as well [2S||. 



B. Dynamic localization 

As we have discussed in Sec. [H] in the strong dynamic 
localization regime the residual absorption is determined 
by dephasing. Using the results of the previous sections 
we can identify three sources of dephasing. 

(i) Escape to the contacts. The quasiparticlc relaxation 
rates for the sequential tunneling and inelastic cotunnel- 
ing are given by Eqs. l(T§l) and (J2SJ. 

(ii) Phonon emission. According to the arguments 
given in the end of Sec. II V Al we can write 



7 P h(T) ~ 6 ( 



(51) 



(iii) Electron- electron collisions. The corresponding 
quasiparticle relaxation rate in a quantum dot was cal- 
culated by Sivan, Imry and Aronov 30] : 




; 0.4 - 



300 



U, (ieV 



(52) 



where Eth is the Thouless energy. The derivation of 
this expression implies the effective continuity of the 



FIG. 3: Electronic temperature in kelvins versus reduced 
gate voltage: (a) only sequential tunneling is taken into ac- 
count, (b) cotunneling is added, and (c) cooling by phonons 
is present as well [2S|. 
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G/Go 



FIG. 4: The function W{G/G ) defined in Eq. for se- 

quential tunneling. 



many-particle spectrum, which imposes a condition T* 3> 
y^ErhS/ \n(ETh/S) 31] . Obviously, for the dynamic lo- 
calization to have any chance to develop, the condition 
7<^(7*)t* <C 1 should be satisfied. 

Suppose for a moment that dephasing is dominated 
by electron-electron collisions, while cooling is dominated 
by the escape to the contacts (later we will analyze the 
conditions for this to be true). One can notice a common 
property of Eqs. ifTBT) . ifTHl) . and for both sequential 
tunneling and cotunneling G/Go and W ut(U)/[(j/S)T 2 ] 
are functions of x = U/ (2T) only. This allows us to write 
a relation 



W 



(j/S)T 2 W(G/G ) 



The energy balance condition takes the form 



Wu 



T 2 

E Th 



Jr 2 W(G/G ) 



(53) 



(54) 



or ( 7 /(5)W(G/G ) = {TjE Th f. Since U and T have 
dropped out, the solution of this equation for G is inde- 
pendent of U, leading to a flat plateau on the Coulomb 
blockade curve G(U) With logarithmic precision 

W(G/G ) ~ G/G (see Fig.gJ, so the level of the plateau 
is G/G^{5h){T,/E Th f. 

Note that the largest possible value of W(G/Go) is 
7r 2 /12 reached at G/Go = 1 (corresponding to U = 0). 
Therefore, the solution exists only if 



Physically, this means that the dot should be sufficiently 
open, so that the cooling is intense enough and the sta- 
tionary temperature is not too high to destroy the local- 
ization. Note that for the observation of the plateau the 
condition 7 <C l/t* is not necessary: even if at U = 
the dynamic localization is absent, as U is increased, the 
dot becomes effectively more closed, so the dephasing by 
escape becomes less efficient. Of course, for the Coulomb 
blockade itself to be present, the condition -f/5 <C 1 
should be satisfied |32 . 

Now let us consider the very top of the peak, [7 = 0. 
Including the dephasing due to both escape and electron- 
electron collisions, we can write the energy balance con- 
dition as 



7 r2+ T 2 
5~ 5 E? 



2 T 2 



(56) 



Th 



Here the left-hand side represents the cooling rate in the 
peak, the first term on the right-hand side comes from 
the dephasing due to escape, and the second term repre- 
sents the contribution from collisions. Due to the condi- 
tion (|55|l the second term is necessarily negligible com- 
pared to the left-hand side, so the only way to satisfy 
the equation is to have T(U = 0) ~ T*. Thus, for the 
dynamic localization to be possible the dephasing in the 
very peak of the Coulomb blockade must be dominated 
by escape. 

Eq. H4tjf> remains valid in the dynamic localization 
regime as well, as it does not depend on the details of 
heating and cooling mechanisms. Thus, one can extract 
the temperature of the stationary state at U = measur- 
ing the curvature of the peak, and study its dependence 
on control parameters: intensity T and coupling to the 
contacts 7. From Eq. I|56|l it is seen that this dependence 
is the strongest when j/S is close to T 2 jE\ h (up to a nu- 
merical coefficient), i. e. when the dynamic localization 
in the peak is about to be destroyed. If we plot T(U = 0) 
versus T (Fig. [SJ , we see that destruction of the dynamic 
localization manifests itself as a crossover from the linear 
dependence T tx T deep in the DL regime (small Y) to 
T oc vT in the Ohmic regime. According to the above- 
said, this crossover can be quite pronounced (like shown 
in the figure) when 



r r 2 1 

ni Th I* 



(57) 



As U is detuned away from the peak, the dot be- 
comes effectively more closed, and simultaneously the 
electronic temperature grows and electron-electron col- 
lisions become more frequent. Thus, the crossover from 
the peak to the plateau occurs where the two mechanisms 
are equally efficient. With the logarithmic precision this 
happens at 



• 

S \E Th 



(55) 



Tl max 



7 E Th \ 

's n j 



(58) 
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FIG. 5: Dependence of the electronic temperature in the peak 
T(U — 0) on the perturbation intensity T for 'y/S = 0.02, 
S = 0.3 /ieV, cu — 3 /ieV, Eth = 100 fieV; in reality the sharp 
angle is replaced by a smooth crossover. 



depending on whether the plateau is in the region of se- 
quential tunneling, j/S <C T*/ ' Ej<h, or of inelastic cotun- 
neling, j/S 3> T^/E^h- The plateau ends when the tem- 
perature of the dot becomes so large that the dynamic lo- 
calization is destroyed by dephasing. Obviously, this hap- 
pens when the horisontal line G/Go = (6 / ln f)(T* / 'Eth) 2 
hits the curve (|48|l or H49JI . which happens at 



U~U n 



E Th 



7 E 



Th 



5 Tl 



(59) 



The resulting shape of the Coulomb blockade peak is 
drawn schematically in Fig. for the Ohmic absorp- 
tion and dynamic localization regimes. The two bound- 
aries H580 and i|59|l give a nonzero range of U (i. e. 
U m in < U m3lX ), if T* <C Exhy/S/T, which can be equiv- 
alently rewritten as 7 e _ e (T») <C l/f», i. e. a necessary 
condition for the dynamic localization itself. 

It is convenient to introduce two dimensionless param- 
eters, corresponding to two experimentally controllable 
parameters T and 7: 



r / w 

~6 \Er^ 



2/3 



\ -2/3 
— ) 



(60) 



The condition 7 e _ e (T*)t* -C 1 becomes / -C 1, the condi- 
tion (|55|l is y 3> J 2 . The top of the peak will correspond 
to DL regime if jt* < 1 or Jj/ < 1. The resulting "phase 
diagram" is shown in Fig. [7| The conditions (|57|l for a 



FIG. 6: A sketch of the Coulomb blockade peak shape in the 
dynamic localization regime without taking into account the 
phonon cooling (solid line): at small U < Umin the dephasing 
is dominated by the electron escape (peak), at larger U - 
by electron-electron collisions (plateau), and finally, at U > 
U m ax the cooling is insufficient, the dynamic localization is 
destroyed, and the dot is in the Ohmic regime. The Ohmic 
curve is also shown for reference by the dashed line. 



pronounced crossover in Fig.JSJcorrespond to crossing the 
parabola y = I 2 in its lower part. 

So far, when analyzing the dynamic localization, we 
did not take phonons into account. Now consider another 
extreme case: both cooling and dephasing are entirely 
due to phonons. Then the energy balance condition in 
the localization regime reads as 



ji4 



Hn6 



= W n 



(61) 



giving T r*> T*. Note that this conclusion is independent 
of the power of temperature in the phonon cooling rate [or 
of a appearing in Eq. (|36[l ] . Obviously, phonons will dom- 
inate if T ph < T*((5/7) 1 / 4 , T ph < y/T*E Th . In this case 
the shape of the peak is given explicitly by Eq. (|15fl . its 
tails - by Eq. I|25[) , and the width corresponds to the elec- 
tronic temperature of the dot. The signature of the dy- 
namic localization effect would be the linear dependence 
of the temperature on the microwave power, in contrast 
to the 1/6 power for the Ohmic absorption case (see the 
previous subsection). The localization regime exists as 
long as 7 P h(T*) < l/£«, or T* 5 < wT p 4 h . The solution for 

the Ohmic regime is T max = (wJ 1 ,^) 1 / 6 , and it is stable 
as long as j p h{T max ) > l/U, which gives T* 5 > wT p \. 
Thus, at a certain intensity such that T* ~ (^T^) 1 / 5 
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FIG. 7: A schematic view of the "phase diagram" in terms 
of the dimensionless intensity and escape rate (7 — y plane) , 
without taking into account cooling and dephasing due to 
phonons. The top of the Coulomb blockade peak corresponds 
to dynamic localization regime only in the region 1; the fiat 
plateau in the tails exists both in regions 1 and 2; in the 
region 3 DL is absent. 

there is a crossover between the localization and Ohmic 
regimes. 

Including all mechanisms, we can note that if the 
electron-phonon interaction is weak enough, T^ h ^> 
E^ h u>/T^ , the phonon cooling plays any role only in the 
Ohmic part of the Coulomb blockade tail. Otherwise, 
phonons start to "eat up" the plateau from the large U 
side The plateau will disappear at T p h ~ \JT*Eth- 
As an illustration, for the intermediate case, we plot the 
Coulomb blockade tail in Fig. [5] in the dynamic local- 
ization regime with and without phonon cooling (lower 
and upper solid curves, respectively) together with the 
corresponding Ohmic curves shown by dashed lines. 



VI. CONCLUSIONS 

We have studied electronic conduction through a quan- 
tum dot in the Coulomb blockade regime under an exter- 
nal periodic perturbation. In contrast to the well-studied 
equilibrium case, the electronic temperature of the dot 
under pumping is different from that of the contacts and 
the substrate. It is determined by the balance between 
heating by the perturbation and cooling due to electron 
exchange with contacts and phonon emission. When the 
cooling is dominated by the former mechanism, its rate 



FIG. 8: Normalized conductance versus reduced gate voltage 
(Coulomb blockade tail): dynamic localization regime with 
(curve a) and without (curve b) phonon cooling taken into ac- 
count, and the same for the purely Ohmic absorption (curves 
c and d) H|. 



depends on the gate voltage, and so does the dot temper- 
ature. As the gate voltage is detuned away from the peak, 
the cooling rate decreases, and the temperature increases. 
As a result, the tails of the Coulomb blockade peak fall 
off less rapidly than in the equilibrium case: instead of 
the usual exponential fall-off for the sequential tunneling, 
under pumping one has a power-law dependence 14811 . 
while for the inelastic cotunncling the equilibrium power 
law is replaced by a weaker one, Eq. I|49[l. At sufficiently 
high temperatures cooling by phonons becomes impor- 
tant, which sets an upper limit for the dot temperature 
(depending on the pumping intensity), which however, 
can be significantly higher than the cryostat tempera- 
ture. 

In the strong dynamic localization regime the heat- 
ing rate is determined by dephasing, as the usual linear 
absorption is blocked by quantum interference. The de- 
phasing can be due to electron-electron collisions, elec- 
tron escape to the contacts, as well as phonon emission. 
The most peculiar situation is realized when the cool- 
ing is due to the contacts, while the dephasing is due 
to electron-electron collisions: in this case the Coulomb 
blockade peak has a flat shoulder, where the conductance 
does not depend on the gate voltage. Such a shape could 
be an experimental signature of the dynamic localization 
effect. 

Finally, we wish to note that conductance measure- 
ments are not necessarily the only possible way to de- 
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tect the dynamic localization. An isolated mesoscopic 
sample can be put into a microwave cavity, and the en- 
ergy absorption rate can be measured as it affects the 
Q- factor of the cavity 34] . In this case the only cool- 
ing mechanism is phonon emission, while the dephasing 
can be due to electron-electron interactions as well. In 
the dynamic localization regime the absorption rate de- 
pends nonlinearly on the ac field intensity: W- ul oc T 3 
if the dephasing is dominated by electron-electron colli- 
sions (T ph ^> ^T*E Th ), or Wi n T 6 if the dephasing is 



dominated by phonons (T pn <C VT^Erh). 
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